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ABSTRACT 


A set of two dimensional finite elements were developed 
to approximate three dimensional coupled plate-acoustic systems. 
Though quite generally applicable, these elements were restricted 
to problems in small enclosure acoustics and sound transmission 
between rooms. 

An experimental apparatus was constructed to substantiate 
the theoretical results. Agreement was found to be good in all ex- 


cept the most adverse cases. 
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CHAPTER I 
INTRODUCTION 


The study of sound, similar to most other forms of 
science, found its origin in a practical "trial and error" sense 
long before the mathematical theories were developed. Experimentation 
on motion by Galileo Galilei [1]* and extensive research on musical 
instruments and listening rooms by Mersenne [2] were grouped under 
the field of "new sciences" until the name acoustics was introduced 
by Sauveur [3]. These early treastises, although mostly non-mathematical, 
defined the field and even made an approach at dividing it into sub- 
fields - geometrical, statistical and wave acoustics. 

Geometrical and statistical are more generally grouped under 
architectural acoustics as they were the first techniques to find en- 
gineering application. The former allows the sound waves only to have 
ray characteristics thereby allowing the path of the sound to be traced 
much as a light ray with its intensity diminishing with distance and 
reflection. Statistical acoustics derives its name from the requirement 
of a "statistically diffuse sound field." Diffuse in the sense that 
the pressure level within the volume can be considered uniform. This 
can usually be approximated in very large enclosures, i.e. auditoriums, 
large halls, etc. where the longest sound wavelength present is less 
than one-tenth the minimum dimension and the multiplicity of sound fre- 
quencies is sufficient to eliminate pressure nodes and appreciable 


pressure maxima. This also allows the wave characteristics of sound 


*Entries in [ ] square brackets refer to the bibliography. 
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to be totally neglected. Instead, sound radiation becomes independent 
of spatial location and the time dependency of a sound is only a 
function of the statistical average of damping surfaces within the 
enclosure. Much of the early work done in architectural acoustics 

[4, 5, 6] took this form of analysis. 

Wave acoustics does not neglect the obtrusive wave charac- 
teristics of sound, increasing many-fold the difficulty of mathematical 
calculations. It is necessary to include reflection, absorption, dif- 
fraction and phase shift at every incidence of the sound waves with a 
boundary, not to mention cancellation and reinforcement of sound waves 
coincident within the enclosure. This approach can be used for any 
size of enclosure but is usually limited to smaller volumes, i.e. small 
rooms, ducting, etc. where the sound wavelengths present are of the 
same order as the enclosure's dimensions. The shear futility of 
analyzing a geometrically complicated structure can be appreciated. 
Major early theoretical contributions in wave acoustics were made by 
Rayleigh [7] and Lamb [8]. Their contributions by no means embody the 
entire field of acoustics, but it has been said, ". . . it behooves 
all of us. . . to read Rayleigh's Science of Sound regularly." [9] 

It might be well mentioned here that between these two 
approaches to acoustics lies a zone that must be approached with an 


extensive "salting" of experience. 


1.1 Review of Related Work 
In 1944, Morse and Bolt published a "state-of-the-art" 


article [10] on room acoustics. A discussion of geometric and wave 
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acoustics was given with special reference to different types of ab- 
sorbing wall materials. These panels were characterized by an acoustic 
impedence which was the complex ratio of the acoustic pressure and the 
air velocity just outside the surface. The real part of this ratio is 
the dissipative term whereas the imaginary part is the reactive term. 
The exact analysis of a regular acoustically hard rectangular room with 
one rigid but absorbing wall was presented but the intricacies in the 
calculations caused the authors to go to more approximate solutions for 
general application. The introduction of more rigid, absorbing walls 
only complicated the analysis. A simplified model that could handle 
different geometries and absorption at the walls would be desirable. 
Principal mode analysis supplies an effective method of 
studying mechanical vibration. The inclusion of damping (hysteretic 
and viscous) cause problems but these are not insurmountable. Bishop 
and Gladwell [11] discuss this method quite fully concluding that 
", . . when damping of the most general type is present, the system 
can vibrate in a principal mode at a natural frequency without there 
being any phase difference between the displacement at various points." 
This would tend to indicate that a study of undamped systems would 
supply an insight into lightly damped systems. The orthogonality of 
the principal mode allows the deformed shape of the structure to be 
represented by a linear combination of the principal mode shapes. 
Though principal mode analysis was discussed in terms of structural 
systems, it could be applied equally as well to acoustics with the 
corresponding "stiffness" and "mass" effects (more correctly kinetic 


and potential effects, respectively). 
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Later, authors became more interested in cases where the 
flexural natural frequencies of the panel were in the same range as 
the acoustical natural frequencies. The previous model of a com- 
plex acoustic impedance for an absorbing panel was no longer suf- 
ficient. Sound transmission between rooms [12] and noise effects on 
flexible structural enclosures [13, 14, 15, 16, 17, 18] can only be 
handled by a coupled plate-acoustic system. The idea of such a coupled 
system, with slight alteration, is also used in the study of liquid 


Sloshing in containers [19]. 


1.2 Aim of Thesis 

Previous finite element work by Craggs [15, 16] produced large 
three dimensional acoustic and plate elements. They were used to study 
transient responses of a plate-acoustic system. Each element required 
relatively large computer core area, therefore even small problems re- 
quired a long execution time which is approximately proportional to the 
Square of the storage space required. By assuming two of the boundaries 
to be parallel and hard, the exact solution can be used for the direction 
perpendicular to these boundaries, i.e. the z-direction. A similar 
approach (sometimes credited to canton ch) was used by Cheung and 
Zienkiewicz [20] in the study of plates and shafts. A reduction of 
element storage area by a factor of 3.0 and a consequent reduction in 
computation time by 7.0 is realized using this assumption. Though the 
elements are now effectively two dimensional, total freedom is retained 
in the remaining coordinate directions as to method of solution. It 
was proposed to develop second order plate and acoustic elements, i.e. 


approximating displacements and pressures and their first derivatives. — 
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Also, instead of a transient study of a coupled plate- 
acoustic system, a principal mode analysis was undertaken. The 
natural frequencies and principal mode shapes of a system can be ob- 
tained by studying the free vibration problem. This entails assuming 
that all time dependent variables obey the same harmonic function, 
thus removing the time dependency of the system. The problem then 


reduces to the solution of an eigenvalue problem. 
{LA] - ALB]} {x} = 0 


The numerical results presented are the natural frequencies 
and principal mode shapes for: 
i) a rigid room with a sloping wall 
ii) two rigid rooms coupled by a flexible panel 
of various widths 
iii) two rigid rooms coupled by an air slit of 


various widths 
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CHAPTER II 
UNCOUPLED ELEMENTS 


The finite element method is basically the approximation 
of the desired solution over a well defined geometric space, i.e. 
cuboid, tetrahedron, etc. Interpolation functions modified by arbi- 
trary parameters make up the general form of the solution. This 
general form is then forced to satisfy the equations throughout the 
well defined geometric space (from now on called an element). Though 
it is not necessary, the interpolation functions are usually picked 
in such a way that the arbitrary parameters take on physical signifi- 
cance at certain points in the element (called nodes). What remains 
is to join several elements into a desired configuration that has 
some degree of continuity between elements and to satisfy the natural 
boundary conditions. 

Convergence (and rate of convergence) to the correct solution 
is of interest. It is extremely difficult, if not impossible, to prove 
convergence to the correct solution for all configurations which can 
be produced by finite elements. The generally accepted approach [19, 
20] is to obtain an approximate solution to a problem that is mathe- 
matically solvable. Successively finer "grids" of elements are used 
to better approximate this solution. Percentage error (between the 
exact and approximate solution) is plotted against the number of ele- 
ments on logarithmic graph paper. The slope of the resulting curve 
is an indication of the rate of convergence. Beyond this, "the proof 
is in the tasting". Results for other desired geometric configur- 


ations obtained from this method must be accepted (with due caution, 
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of course) and if possible checked out by experimentation. 


2.1 Element Formulation 

It has been found more convenient to work with functionals 
than with the actual equations of motion and the boundary conditions 
because certain symmetries arise in the functionals that can be taken 
advantage of. If a functional can be found whose stationary condition 
Satisfies the equation of motion and the natural boundary conditions, 
then the functional can be used just as effectively in the solution of 
the problem [21]. In fact, the functionals required for the uncoupled 
panel and acoustic subsystems are the Hamiltonians or action integrals 


J of the subsystems 


ty 


ae isa lee | dt 
oT 
where T is the kinetic energy and U is the potential energy. The 
governing equation and boundary conditions are obtained by setting the 
first variation to zero 
6J = 0 

The first variation is obtained by replacing the required 
parameters by themselves plus a small arbitrary increment (i.e. p is 
replaced by p + 6p) then the resulting equation J + 6J and the original 
equation J are subtracted to give the first variation. 

Appendix A offers a relatively detailed formulation of the 
acoustic finite element. From the basic acoustic and energy equations 
(Appendices Al and A2, respectively) through to the finite element or 
matrix formulation of the acoustic functional (Appendix A4), there 


was no assumption made as to the form of time dependency. This would 
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allow the study of transient acoustic systems; however, since the 
steady-state response of an undamped system is desired, the time 


function was written as 
sin «© t + 6) 


where w is the cyclic frequency and 6 is a phase angle defined by 
initial conditions. Since all variables in the system obey this time 


function (Al16)* of Appendix A4 can be written 
(sale oc (Ped? 1p 1. = 0 
e ee ne 


which is a standard form of the eigenvalue problem with pve corres- 
ponding to the elemental nodal pressure coordinates. For a system com- 
posed of several elements, the ‘e' subscript can be dropped. 

The overall [S] matrix can be derived by noting that the energy 
of the total system is a linear summation of its component energies [8, 
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with a similar formula for [P]. Leading to 
{[S] - #[P]} {p,} = 0 (2.1) 


where {p,} is the overall pressure coordinate. 


*Entries in ( ) refer to equations. 
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Figure 2.1 (a) RECTANGULAR ACOUSTIC ELEMENT 


Figure 2.1 (b) TRIANGULAR ACOUSTIC ELEMENT 


Figure 2.1 (c) LINE PLATE ELEMENT 


Figure 2.1 ELEMENT SHAPES AND NOTATION 
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2.2. Rectangular Acoustic Element 
Figure 2.1(a) gives the element shape and dimension in the 
x-y plane. The encircled numbers are the nodal numbers and the nodal 


values of p, Sp, Sp are p., Sp., Sp., i = 1, 2, 3, 4 which are functions 
dies HB hy Aa ae 


of time. Since this is a "pseudo" three dimensional approach, there 
Should be an identical x-y plane at a distance d behind the plane shown. 
These planes were previously defined as rigid and parallel, allowing the 
exact solution for the z-direction to be found. This solution, a single 
Fourier cosine series, can easily be used to approximate many mode shapes 
in the z-direction (as is shown in Appendix B1), although only the nth 
term of the series will be used in the following calculations to simplify 
the notation. In the x-y plane, pressure and pressure variation can be 
approximated by the product of two sets of Hermitian interpolation poly- 
nomials f(x*), f(y*), where x* = x/a and y* = y/b. These polynomials, 
their characteristics and their integrals are given in Appendix B2. The 
desirable characteristics of these interpolation polynomials allow Pas 
the pressure distribution throughout the element for the nth node in the 


z-direction, to be expressed as 
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where the dimensions a, b are a result of the chain-rule, i.e. 


ope ox* or 6p = a Op. A simpler short-hand notation will be 
2 oe on 


introduced here 


mi 


Pe = {f(x*) f(y*) {p Je COs a8 


where the superscript "T" indicates the transpose of the vector and 


Fy (x*) Fy (y*) bp, 
afa(x*) £,(y*) ae 
{f(x*) f(y*)} = : and bis : 
bf4(x*) fy (y*) So, 


It would be best to note here that 
8, = 1 te) f(y*)} fp}, cos ntz 
pa n e eee, 
X a d 
Pe = 1 {f(x*) f' (yt) } fp Peecose nilz 
— n e ae 
y, b d 
Sp = <n {F(x*) f(y*) fp Jasin nz 
Z d d 
The prime "'" denotes differentiation with respect to the variables in 
round brackets. 
Substituting these formulas in the energy equation as done 
in Appendix A4, the acoustic stiffness [Sep and pressure [Pode matrices 


can be written 
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The integrals with respect to z are defined in Appendix Bl, 
and the integrals with respect to x* and y* are defined in Appendix B2. 

The rectangular acoustic element was programmed for the IBM 
360/67 computer and with the aid of the IBM Scientific Subroutine Pack- 
age, the eigenvalue problem, as shown by (2.1) was performed. A cubical 
room of unit side and rigid boundaries was approximated as a check for 


convergence. This system has an exact set of principal mode solutions 


p = p* cos f~mx * cos mmy * Cos nnz 


where the pressure p* is a periodic function of time and 2, m, n are the 
number of half-cosine waves between the boundartes in the x, y, z dir- 
ections, respectively. Grid refinement was done only in the x-direction. 
Percentage error in the eigenvalue (the square of the natural frequency) 


versus the number of elements is plotted for various (£, m, 0) mode 
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Figure 2.2 CONVERGENCE OF RECTANGULAR ACOUSTIC ELEMENT 


13 


2inemela te vedmui 


f 7 7 » 
THIM313 ITZUOOA AAJUOMATIIA 70 3OMIOAIVMOD S.S 


14 


The convergence of the results is at the rate of ine, 

(n, being the number of elements) for (2, 0, 0) modes because of the 
x-direction grid refinement. The two (f, 1, 0) modes plotted, which 
had no y-direction grid refinement, show no improvement with increase 
in element number. Note that even one element made a good initial 
approximation. Although not plotted, (f, m, n) modes would give ex- 
actly the same type of curves as (f~, m, 0) curves because the exact 
solution between rigid boundaries is used in the formulation of the 


elements for the z-direction. 


2.3 Triangular Acoustic Element 


Because of the change in geometry, the method of approxi- 
mating the pressure Po within the element was done by a cubic poly- 


nomial in x and y, and the exact solution in the z-direction. 


3 3-m Lom 
Po =( am X Y ) cos nz 
m=0 p=0 d 


where the Bam were a set of arbitrary constants. 

A complete cubic polynomial has ten arbitrary parameters, but 
to approximate the pressure and the pressure variation in x and y- 
directions for three nodes required only nine parameters. There are 
several methods to remove the extra parameter; two were tried but the 
one presented in Appendix Cl was rejected, in this case because the 
added computational time required did not improve (and in some cases 
worsened) the element's accuracy. The method chosen was to simply 
equate the Ao] and aio coefficients, so the elemental pressure Po 


was written as 
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The problems with this formulation will be discussed later. 
A third possibility would have been to define an internal node where 
pressure was the only parameter, but this would cause problems with 
the storage and manipulation of the resulting matrices. 


Using a notation similar to Section 2.1, 


Daa {g(x,y*) }! {atcos nmz 
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{g(x,y)} = ]y fand {a} = Jay 
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Considering the element and notation of Fig. 2.1(b), it is 
now necessary to relate the arbitrary constant {a} to the nodal charac- 


teristics Pas Sp. 5 Sp. > i= 1, 2, 3. This is simply done by writing down 
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where 


Pate 5 oP 


Sp. 
Sys 
and [A] is a matrix of the nodal values of {g(x,y)} and its directional 
derivatives {sa and {sah » respectively. Then it can be written 
Ox éy 


fe) I ate 


If the formulation was carried to completion from this point, 
the resulting integrals over the triangle would be arduous and could 


only be handled numerically. It is more convenient to set up a localized 


x' - y' coordinate system, for the triangle, where x' = e€ + b and 
a+b 
y' = n/h perform the required integration (shown in Appendix B3) then 


"rotate" the results into the x-y coordinate system (shown in Appendix 
B4). The formulation for the x' - y' coordinate system is identical 

to what has been presented for the x - y coordinate system, but gives 
some simplifications in [A] and Ale as shown in Appendix B5. The 
formulation can be carried on in terms of x' and y'. 


The elemental pressure can be written as 


pe = t9(x'sy') CAI”! fp}, 


fasofsosath ett bas {(y.x)p} to esulsv isbon a. 


eat 


xfbnaqqA nt nwode) moctave apiece 2 Xx sft ont estas 


iso0i 6 qu dsz ot shetasvN0D ovo eft $1 eT soiyoiue bef baad 


oAT .24 xtbnoqgA ni awode 2s ,/ [A] bne CA] AF ane 


ods 7 
natitw sd asd Jt naHT ah e " ‘as 
es 


or 
roq efds movt noftslamoo ot batavs> 26w ~thenual mil 


bfuos bre evoubys ad Bluow slpnstyd sat 4evo ctevpstnt tat 


- 
ba6 d+ 2 = 'x sysdw -Slonsryd sdt vo? emoateye 9 


ois (£8 xtbaeqqA nf nworiz) no fSeyposn? bervtupsy ade t 


lsorinebS et mateve stsnfbxo005 'Y - 'x ome OF nottety } 


Pin oa 


. . ; 5 as 
fe Jud .mosevye atentbioos y - x Bat vor bsinezeyg « 


ta 

‘\ bis 'x to ems at ao Series 6d neo aOre 
, | . ; es 

26 NSTIiaw sd nso syuzeong T stnoma fe ent 


3f4a) | CAD. 74) We xB} = 39 


Note that the variables in round brackets are omitted in future to 
Simplify the formulation. The acoustical stiffness [S,'], and pres- 
sure [Po]; matrices are derived by substituting into the energy 


equations (Al14). 
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The integrals with respect to x' and y' are given in 
Appendix B3. The integrals with respect to z are given in Appendix 
Bl. The rotation of the nodal velocity potential {on}, from the 


x' - y' coordinate system, is shown in Appendix B4. Since energy is 


independent of coordinate transformation, it can be written 
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where [6], is the acoustic rotational matrix defined in Appendix B4. 

The triangular element acoustic stiffness and pressure 
matrices derived above were programmed and the convergence tests were 
run on the cubical room used in Section 2.2. The orientation of the 
elements and the percentage error versus number of elements plot, are 
Shown in Fig. 2. 3. The rate of convergence was slightly less than 
1/2 in the direction of grid refinement with a poorer initial approxi- 
mation than the rectangular elements, but overall still gave good 
results. It was interesting to note that by changing the pattern of 
subdividing the cuboid, even for a one-directional grid refinement, 
the results (and the resulting percentage errors) were changed, but 
this could be a study in itself, so will not be gone into here. 

By studying the formulation of the exact solution and by 
comparing the finite element equations (2.3) and (2.5) it is possible 


to arrive at a natural frequency equation for the (m, n, q) mode in 
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terms of the (m, n, 0) mode. 


2.8) 


where Ung is the natural frequency for the (m, n, q) mode. 


c is the velocity of sound 


d is the depth of the room 


2.4 Beam Plate Element 

A similar detailed formulation to the one given in Appendix 
A could be undertaken for flat plates, but would serve no useful pur- 
pose, so will be forgone in favour of writing only the energy equations 
according to S. Timoshenko [23] for cartesian coordinates. It is, 
however, initially necessary to discuss the orientation of the plate. 
To use Kantorovich's approach the plate must be placed between the two 
rigid planes perpendicular to the z-direction, but the other direction 
must be general. It will be seen that the same rotational concept 
used in the triangular acoustic element can be used here. Therefore, 
the plate will initially take an x-z plane orientation for the formu- 


lation and later be rotated to the x' - z'plane by a rotational matrix 
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where Wo is the time derivative of the transverse deflection of the 
plate and dxdz is an infinitesimal surface area on the plate in the 
X-z plane. 


It is also convenient to introduce the differentiation 
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2) is the Plate Modulus 


ff Seren ein 
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E is Young's Modulus 
v is Poisson's Ratio 
Again, it is necessary to show that the equation of motion 
and the natural boundary conditions can be derived from the action 


integral by variational calculus. 
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The rest of the terms in (2.11) can be integrated by parts 
twice, with due care in the evaluation of the integral. The third 


term in (2.11) can be written as 
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where the first and second terms are evaluated at the front and rear 
boundary (z=0 and z=d, respectively) before they are integrated with 
respect to x. The last term in (2.11) is divided into two halfs and 
integrated twice by parts; one half is integrated first with respect 


to x, then with respect to z; the other half is integrated in the 
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reverse order. Noting again that SW) and SW 2x are arbitrary and 


small and SW, = 0 at ty and tos the following equations result. 
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These sets of equations show the power of the functional 
approach. Besides (2.14) being the equation of motion for a free 


vibrating plate, (2.12) and (2.13) embody all the natural boundary 
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conditions on a rectangular plate. According to Timoshenko [24] 


Mxd. OF (Wy xx ' Wy 222) 
Le = - D (1 - v) Wo 2x2 
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where Ne is the bending moment (in the x-direction on the plate edge 
perpendicular to the x-direction), ee is the twisting moment and OF 
is the shear on edges perpendicular to the x-direction. 


Using these relations in (2.12) 
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For example, a pinned edge on x = O is defined by ow, = 0 and 
M.. = 0. This allows for non-zero Qy and Moy: All other boundary con- 
ditions [24] can be derived from these equations. 
Although (2.14) has no assumption as to the time function, 
it will be considered as harmonic (the time function) and real, i.e. 
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For simplicity the exact solution in the z-direction will 
be taken as a single Fourier sine series, i.e. pinned boundary con- 
ditions at z = 0 and z = d. Only one term (the kth) of the exact 
solution is to be used to simplify notation, although again a multi- 
term approximation of a mode shape could be handled. 

If the notation used in Fig. 2.1 (c) is adopted (noting that 
6 = 0 for initial consideration), we can write the deflection of plate 
w., in terms of a product of a set of Hermitian interpolation poly- 


p 
nomials {f (x*)} and the exact solution where x* = e/fy and O<e<f,. 


Wo = {Wy fF, (x*) + Ly F4(x*) 


Wo F4(x*) + is 2s Fy (x*)} sin ee (2.16) 


where w., SwW., 1 = 1, 2 - are the nodal characteristics which may vary 
OX 
with time 


F(x), i= 1, 4 - are the Hermitian polynomials 
definition Appendix B2. 


This can be written in short hand notation as in Sec. 2.2 
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The elemental mass [MJ and stiffness [KJ can be defined 


by the equations 
T, = 1 tw}! [MJ] tH} 
e 0 n e n 
ul=14w! EK. 4a} 
ely n e n 
By substituting into the respective energy equations (2.9) and (2.10), 
[K,] = Dg ; CLF! (x*)} CFE! (xe) Hy et CF(x*)} (FO) }) fk : 
Poe Lp ca 
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0 d 
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where the third term has been arranged to retain the symmetry of the 
stiffness matrix. 
aa 
[MJ = eph Ly ; CF(x*)} C(x)! dx* +f sin’ knz dz (2.18) 
Bu 0 d 

The integrals with respect to x* can be found in Appendix 
B2. The integrals with respect to z can be found in Appendix Bl. 

This element can now be rotated to any position, letting 


1* denote the dimension along the sloped plate as shown in Fig. 2.1(c). 
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By direct comparison with Appendix B4, 
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The plate rotational matrix [6] can be written 


Sin 0% cosré 0 


O 
0 0 ] 
ca, = 
: sin 0 COS 10 0 
O 
0 0 ] 


These plate element stiffness and mass matrices were pro- 
grammed for the IBM 360/67 and convergence tests run on harmonically 


vibrating plates of unit side with pinned boundary conditions on all 


edges. This plate has an exact principal mode solution 


Wo = Wye sin rmx sin ktz 
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where r and k are the number of half-sine waves in the x and z- 


direction, respectively. 
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au, ho and D were picked so the coefficient in front of 


Cal (r° H Ko), 


p 
the natural frequency equation was unity, i.e. w 


There is no similar formula to (2.8) that can be easily 
derived for a plate system so the plots of percentage error versus 
number of elements are given for both k = 0 and k = 1 (Fig. 2.4 and 
Fig. 2.5, respectively). The initial approximation in some modes are 
about 60 per cent in error in the square of natural frequency 
(8 per cent in the natural frequency). This is not good, but for- 
tunately a slightly higher number of elements gives accuracy of 
2 per cent in the lower natural modes. These lower modes will be used 
in the study of sound transmission. It may also be noted that the re- 
sults of the theory are integers close to unity for which the error 
associated with truncation could be significant. Later results tend 
to bear this up. 

The plate element, as developed, was effectively a beam. 
Independently they can be used to study the vibration of rectangular 
plates with two parallel pinned boundaries while the other two 
boundaries can have any combination of free, fixed or pinned boundary 
conditions. Warburton [25] has produced a very comprehensive analytic 
paper on such plates. Melosh, Papenfuss, Clough, Adini and Tocher 
have produced [21] similar plate elements of more general application 


so the use of the plate elements will be left for the coupled system. 


2.5 Application of Acoustic Elements 
2.5.1 Odd Shaped Room 


A problem was required to test a combination of triangular 
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and rectangular acoustic elements. The Mechanical Engineering 
Building on the University of Alberta campus, which was under 
construction at the time, had a reverberation room which offered 
an excellent model. Being placed directly under a sloped roof section, 
required that one half of the room have a sloped ceiling. The exact 
profile of the room is shown in Fig. 2.6 (c), its depth being 
22.058. 

One problem, exemplified by this room, is what boundary con- 
ditions exist at a corner such as shown in Fig. 2.6 (a). Any node 


along the sloped edge has the boundary condition 


Op i="="Ctn Weep “or ~op = =“tan op 
6X éy oy OX 


while along an upright edge 


Right at the corner node, there is a discontinuity, as there 
would be in a mathematical analysis (unless certain transforms were 
undergone). Using finite element methods there are three logical 
possible choices of how to approximate this boundary condition as 
shown in Fig. 2.6 (b). The first two are not correct, but on element 
refinement they will converge to nearly the correct result. This 
occurs by virtue of the fact that as the elements get smaller the 


effect of the elements containing the incorrect nodal value diminishes 
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Figure 2.6 (b) APPROXIMATING BOUNDARY CONDITIONS 
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Scale: 1 in.= 6 ft. ————— coarse grid 


Depth: 22 ft. ———— 1st grid refinement 


Figure 2.6 (c) FINITE ELEMENT GRIDS 
Figure 2.6 SLOPED CEILING REVERBERATION CHAMBER 


Table 2.1 CONVERGENCE OF FOUR LOWEST EIGENVALUES 
FOR SLOPED CEILING REVERBERATION CHAMBER 
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in terms of the overall acoustical kinetic and strain matrices. 
Method 3 could converge to the exact result, but it entails larger 
partitioning operations on the overall matrices (described in 
Appendix C2) and is therefore more prone to numerical divergence. 
If looked at from a potential energy viewpoint, 1 and 2 are respectively 
Slightly weaker and slightly sérohoer than the correct boundary con- 
dition 3 giving respectively lower and higher eigenvalues than 3 for 
any particular grid configuration. Since the approximate solution 
converges on the exact solution from above, it was deemed interesting 
to see the effect on convergence of using boundary condition 1. 
Because of its "relaxing" effect on the overall energy matrices, it 
is reasonable to assume that for a relatively coarse grid pattern, 
boundary condition 1 will give eigenvalues closer to the exact solution 
than a finer grid using the other boundary conditions. The concept of 
"flexibility" of a system is quite often used to describe deviation 
from classical results. A more flexible system has less ability to 
store strain energy and thereby results in a lower natural frequency. 
Grid patterns for the coarse grid and refined grid are 
shown superimposed on one another in Fig. 2.6 (c). There were other 
grids used including a finer one that divided the sloped surface into 
three equal sections, but these would only complicate the figure 
given. The effect of boundary condition 1 is most exemplified for 
the coarse grid where the entire slope edge has the "relaxed" boundary 
condition. For computational convenience c, the velocity of sound, 
was taken as unity. 


By close examination of each column of Table 2.1 it is easily 
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seen that the results for the slope ceiling room are very close to 
the exact solution for a hard rectangular room that would bound the 
reverberation room. This is logical since the reverberation room 
does not deviate much from the classical rectangular form. Column 
one shows the (0, 0, 1) mode, which is independent of the actual 
profile of the room, and should therefore converge exactly to the 
classical result. In accordance with the energy approximation method, 
the results actually start higher than the correct value, drop below, 
then start upwards again. This can be explained by the "flexibility" 
of the structure; the more "flexible" the structure, the lower the 
eigenvalues. Initially the grid was not flexible enough to approxi- 
mate the mode shape, but after the first refinement the structure 
became overly flexible because of b.c. 1. Ensuing refinement would 
lead to closer convergence because the system has sufficient flexi- 
bility to closely approximate the mode shape, while each further sub- 
division of the sloped ceiling lessens the effect of b.c. 1 (until 
One enters the region where numerical errors become significant). 
Similar analysis can be done on the remaining three columns, noting 
that the sloped roof has almost no effect on the (x, 0, 1) modes but 
is quite important in the (0, y, 1) modes. 

Figures 2.7 through 2.9 show the principal mode shapes of 
three rooms of different dimension but with the same basic shape. The 
(a) figures are the actual profiles for the reverberation chamber, 
whereas the (b) and (c) figures show the changes due to substantially 
increasing the effect of the sloped wall on the x-direction and 


x-y directions, respectively. The points of no pressure variation 
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within the room are of importance when making sound measurements. 

For a hard rectangular room the zero-pressure nodal line occurs at 
positions that are simple fractions of the overall dimension, i.e. 
1/2, 1/3, etc. Figure 2.7 shows how the zero-pressure node line 
deviates from the classical position for a rectangular room. The 
narrowing of the enclosure at one end changes the volume distribution 
as well as causing localized compression and rarefaction of the sound 
wave near the roof. A wave travelling in the x-direction, is given 
an unbalanced y-component of force, tending to bend it into a path 
where uniform compression (or rarefaction) occurs on both ceiling and 
floor. In Fig. 2.7 the reflection off the narrow end wall, floor and 
ceiling, at non-zero angles of incidence, causes "focusing" of the 
pressure waves in the lower right-hand corner of the room, giving 
maxima in pressure levels in this region. 

Figure 2.8 is actually very similar to Fig. 2.7 except there 
are now two node lines, i.e. two maxima and a minima, that are sub- 
stantially affected by the sloped roof. It is interesting that the 
maxima along the full height wall in Fig. 2.8 moves from the lower 
left to upper left corner as the slope angle of the roof, n, becomes 
more than 45°. Figure 2.8 (c) has a higher eigenvalue than Fig. 2.9 
(c) so it could be said that the diagram ordering is incorrect, but it 


is of greater interest to show the development of a mode shape. 
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Figure 2.7 FIRST MODE - ROOMS WITH SLOPED CEILING 
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CHAPTER IT] 
COUPLED ELEMENTS 


The elements produced thus far have their independent 
applications, but have not yet been intercoupled, i.e. a displacement 


of the plate must cause a pressure change in the fluid and vice versa. 


3.1 The Coupling Matrix 
Former work by Warburton [26] and Pretlove [27] approached 


coupling by using the exact acoustic equation with a normal mode 
representation of the plate. Zimmerman and Gladwell [28] did a very 
complete treatment of coupled plate-acoustic systems, but encountered 
difficulties because they formulated the acoustic and plated elements 
entirely in terms of displacement or force quantities which, if nothing 
else, cause conceptualization problems. Using the pressure acoustic 
formulation and the displacement plate formulation described by Mason 
[29, 30] and the coupling concept advanced by Craggs [31], the ensuing 
formulation is easily conceptualized and relatively efficient. 

To satisfy equilibrium, the forces due to the flexed panel 
and the pressure distribution over the surface area of the plate must 
be conserved. Terms whose first variation satisfies the boundary 
condition at the interface must be included in the respective 
functionals. This boundary condition is simply that the normal veloc- 
ity of the fluid must be the same as the transverse velocity of the 


plate over the entire plate surface 
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where u + nN is the normal velocity component of the fluid. By 
differentiating with respect to time and using equations (A8) and 


(All) 


grad P) ° n= -0 We = -o grad $ 7 (3) <1;) 
where p is the acoustic density and o, of course, is the velocity 
potential. Recalling the form of the acoustic functional, and adding 
to it a term whose first variation with respect to o satisfies the 
dynamic boundary condition 


is 
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J= sf pe [SSF {grad > + grad ¢ 
ty 2 V 
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eh 
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- p dt 
aS Wo? ds J 
(fe 
where So is the interface surface. This last term can be physically 


interpreted as the pressure displacement work, Bo» at the coupling 


interface 


2s na e > We ds. = ae Py Wp dS, (BieZ) 


It must have a stationary value and satisfy the dynamic boundary con- 
ditions. This can easily be seen by following the procedure of 


Appendix A3 including the extra term. 
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The energy added to (or extracted from) the acoustic element 
in the form of work, must come from (or go to) the adjacent plate 
element. Therefore, the work term must occur in the energy equations 
for the plate, but with the opposite sign. External forces, other 
than the acoustic pressure, can act on the plate elements but they 
will be excluded from this study because it is the free vibration pro- 
blem which is of interest. 

To form the coupling matrix proper, 0, it is now necessary 


to recall the approximations for pressure and displacements used in 


Chapter II. 
Bass tne) f(y*)}! {p,} cos niz (373) 
d 
where 
Ker hia 
ye =: y¥/D 


for the rectangular acoustic elements 


Pp, & tg(x'.y')} {p,} cos nnz (3.4) 
d 
where 
X=) (eve D) 
(a + b 
Va oe eis) 
h 


for the triangular acoustic element. Note the coordinate system 


indicated here is the localized system because this will be the system 
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the plate will be associated with initially. 


Ww. = {#( pry}! {w.} sin knz 
P ; “d 


where 


Depending on the element to be coupled to the plate element, 
one of the two pressure approximations will be used. The pressure 
formulation for the rectangular acoustic elements will be used in the 
following excursion because the experiments and examples to follow are 
based on this element. 

B. = ap Py dS. 

For this example, assume that the line plate element is 

coupled to edge 1 - 3 of Fig. 2.1 (a) where y and y* are identically 


zero. Also, g* = x* and € = x, respectively. So 


. 
B, = {pq} LS {f(x*) #(0)} {F(x#)}! ax* J 
0 p 
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ff cos mnz sin kiz daz) tw} 
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where 0 is the integral within the square brackets. 
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The integrals with respect to x* can be found in Appendix 
B2, while the integrals with respect to z are in Appendix B1. 

Careful examination of Appendix Bl shows that the integrals 
are taken over one-half harmonic interval. This alters the ortho- 
gonality of the sine - cosine product from the standard full interval 
form of Fourier [32], whereas the sine - sine and the cosine - cosine 
product integrals are unchanged. It is obvious then when n + q is an 
even sum of integers, there is no coupling integral produced by this 
method. This would tend to indicate that in the z-direction an odd 
order acoustic mode would be driving an even order plate mode or vice 
versa. This is readily explained in such a coupling situation because 
for n + q = even integer, the net p - v work is zero giving an un- 
excitable configuration. 

The integrals with respect to x* (y* = 0) were for coupling 
for side 1 - 3 of the rectangular element. These integrals will alter 
when another edge of the element is used, i.e. y* = 1, x* = Oor 
x* = ] but all can be handled by the aforementioned Appendices. The 
coupling formulation for a triangular acoustic element is the same; 
substitute (3.4) and (3.5) in equation (3.2) and multiply out the 
entries of matrices for a surface. Because of the random orientation 
in the x - y plane of any edge, it is exceedingly difficult to do the 
ensuing line integrals. It is initially chosen that the surface to 
be coupled, is co-incident with the localized x-axis, as shown in 
Fig. Bl. The integration can then be performed with relative ease as 
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These integrals can be done by hand (some of which are shown 
in the latter parts of Appendix B2). The resulting localized coupling 
matrix Orn] must be rotated by pre-multiplying by an acoustic rotational 
matrix [9], and post-multiplying by a plate rotational matrix [6], 


before it can be placed in the overall coupling matrix. 
[o-prd = £0)! Leotard £0] 
TRI a SORT p 


This completes the formulation of the coupling matrix. 
Including the present discussion, the coupled energy functionals J, and 
Jn? for the overall system can be written from equations in Appendix 


A4 and Chapter 2, Section 4. 
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By taking the first variation of both functionals (the 
first with respect to > and the second with respect to Wo) and 
equating them to zero, the coupled equations of motion can be arrived 
at. To give more physical significance to bd, = 0, this functional 
variation should be differentiated with respect to time and the pres- 


sure formulation substituted. Then, 


i 
S) 


[P] {pt tyis} ap} - [oJ {Wi} (RS) 


[M] {i,} + [kK] {w,} + [o]' tp} = 


I 
(>) 


(3 7) 


These equations of motion are for a coupled plate-acoustic system with 
no external force {F} acting on the plate elements. These external 
forces could be included by simply placing them in nodal form on the 
right hand side of (3.7). {p} and {w} can take any time function desired, 
SO random vibration studies could easily be undertaken. A finite dif- 
ference approach could be used for stepping through time. 

Since this is a steady-state study, this discussion will not be 
pursued. Equations (3.6) and (3.7) can be written in a more compact 


matrix form 


- 
OQ aess -O P p 
where the steady-state time function has been considered to have a fre- 


quency w. The square brackets around the entries have been dropped for 
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convenience. This equation can be used in eigenvalue solutions but 
the matrices are unsymmetric, therefore, computation becomes 
tedious even by numerical methods. Irons [33] advanced an approach 
that makes these matrices symmetric. By partitioning the matrices 
and performing some basic matrix operations giving 

] Te-l 


Ke 80 . M+eS 0 0S P 


= 
0 Pp ps7! ps~!p 


The result is symmetric due to the symmetry of M, K, P and S. This 
approach is not exactly suitable to this situation, as the M, K 
matrices are of the order of one-tenth the size of the S, P matrices. 
To do the inversion of the S matrix would therefore be exceedingly 
time consuming, not to mention inaccurate (because of the small values 
of the entries in the matrix). Instead, the M matrix was inverted 

and Iron's basic procedure was followed giving 


KM Ik omnia! 21m ol}(w 


yO lie best Tak te 0 Pillip 


6M 
which again is symmetric and can be used for the eigenvalue problem 
more simply than its unsymmetric counterpart. 
Initial checks were done on a simple coupled system; a 
cubical room with one flexible boundary (as shown in Fig. 3.1). To 
simplify things even further the flexible boundary was made almost 


rigid (i.e. plate eigenvalues are well above acoustic eigenvalues) 
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so as to deviate only slightly from the exact rigid room case. The 
computed results of the rigid boundaried room and the room with one 
slightly flexible boundary, are shown in Table 3.1, along with the 
exact results for a "rigid" room. The acoustic element is the same 
grid for both computed cases. 

There are two noteworthy items in the results. First, when 
searching for an exact zero eigenvalue, the IBM eigenvalue - eigen- 
vector routine NROOT sometimes gives a negative result. This is 
unrealistic for the positive definite problems for which this numerical 
method is designed. These negative roots must be considered as a 
measure of the numerical error in the calculations. Secondly, in all 
cases, the eigenvalues are lowered by the introduction of the slightly 
flexible wall. Considering only the acoustic equation of motion, this 
would indicate a slight increase in the potential energy (associated 


with [P] matrix) with no effective change in the kinetic energy. 


a? = cw) [s](y) 
{py} [P]{y} 
where {wv} is a principal mode shape. 

From this point on it can be said that there are no further 
independent plane modes, i.e. (x, 0, 0), (0, y, 0), (0, 0, z) modes. 
This is indicated by the reduction of eigenvalues (square of natural 
frequency) for cases in which the plane wave is parallel to the flexi- 
ble boundary. It can also be shown by inspection of the plate eigen- 
vectors, which show deflected plate positions for the above mentioned 


case. The independent plane wave notation will still be used in 
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following discussions, but only for descriptive purposes. 


3.2 Sound Transmission Between Panel Coupled Rooms 


The culmination of all foregoing work would be applicable 
in the study of plate-acoustic systems such as general floor areas 
(with a rigid floor and ceiling) enclosed by or containing flexible 
floor-to-ceiling panels. This is reasonably broad, and therefore 
the simple case of direct sound transmission between two rectangular 
rooms through a variable width floor-to-ceiling panel was carefully 
analyzed and checked. The widths of the rooms were left unchanged 
but the length of one room was altered. An experimental apparatus 
was constructed to approximate the theory. All "rigid" surfaces 
were made of 3/4 inch fir plywood capped on the inside surfaces with 
1/8 inch masonite. A 4 inch speaker was attached to the source room, 
acoustically exciting the interior of the room through an orifice. 
Precision oscillators, frequency counters and pressure amplitude 
monitors were used to measure the standing wave patterns within the 
receiving and source rooms. The source room had fixed dimensions of 
7 inches x 19 inches x 23 1/2 inches depth, while the receiving room 
had similar cross-sectional dimensions but of variable depth. The 
components of the apparatus are shown in Fig. 3.2 (a). Figure 3.2 (b) 
shows the "driver" for the source room as well as the basic hard wall 
construction of the rooms. Figure 3.3 shows the basic layout of 
equipment. 

Detail "A" in Fig. 3.3 shows the mounting frame for the flexi- 
ble plate. To approximate a zero deflection boundary condition, a length 


of 1/2 inch x 1/2 inch aluminum bar stock was accurately machined to 
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Figure 3.2 EXPERIMENTAL APPARATUS 
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accept all edges of a 1/8 inch aluminum plate to a depth of 

1/16 inch. To assure no air leaks, the bar was carefully fitted 
to the rigid plywood wall and silicon sealant was used at all con- 
nections. Though not shown, the upper and lower bar stock rails 
were also fitted to the upper and lower panel with sealant. A 
well defined fixed or pinned boundary condition is not feasible by 
this construction. It was, therefore, necessary to give both the 
fixed and pinned analytical result. 

Detail "B" is the edge connection between the rooms. Firstly, 
the rigid interior wall (or aluminum frame, depending on the flexible 
panel's width) was completely sealed to stop air leakage from one 
enclosure to another. The sealant and the saw cut isolated either room 
from structure-borne vibrations in the exterior walls of the other 
room. This eliminated a large amount of flanking transmission through 
the exterior boundaries which are not ideally hard. The cord binding, 
used to align the rooms and accurately seal the boundaries, was suf- 
ficiently limp to transmit a minimum amount of vibration. As will be 
seen in later results, this isolation was ineffective, in a few cases, 
but on the whole was very satisfactory. 

On a more general scale, the entire apparatus was isolated 
from the laboratory tables by 3 inch foam mats. All boundaries where 
leaks could occur were either caulked and glued or temporarily caulked 
with silicon sealant. The holes for the pressure probe were sealed 
with brass plugs when not being used directly in tests. 

It is of note that the various plates were chosen so that 


their fundamental modes were between 300 ie and 500 Ho which is well 
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within the range of the first few acoustic modes of either room, 
i.e. the plane standing wave parallel to the plate surface in the 
y-direction was 355 Ho. 

3.2.1 Theoretical Considerations 


To obtain accurate values of po, (plate density) and c 


p 
(velocity of sound), careful measurements were made in the test labor- 


atory. They are, respectively 


p, = -003217 slug/in> (for aluminum) 


c = 13,560 in/sec 


Other constants were accepted at their textbook values 


6 


1.376 x 107° slug/in® 


1S 
it 


6 


E.= 10.3, x 10 1bf/in? (for aluminum) 


vy = 1/3 (for aluminum) 


These measurements proved to be valuable because the errors 
associated with published constants are of the same order as the errors 
due to the finite element approximation. 

Se sc RESULTS 


The description of results is broken into two parts: 
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Se) 
i) Effect of changing panel area and stiffness 
on unchanged room volumes. 
ii) Effect of changing length of receiving room 


on unchanged source room and panel. 


Figure 3.4 through Fig. 3.8 are the five lowest non-zero 
modes for identical source and receiving rooms. Each of the diagrams 
in any given figure shows a different panel stiffness and surface area. 
These panels represent 100 per cent, 73.7 per cent and 42.5 per cent 
of the coupling wall, respectively. Because the rooms are symmetric 
about the panel, the principal mode shapes are symmetric, therefore, 
only one half of the mode shape will be plotted. The positive and nega- 
tive signs indicate the relative pressures and rarefactions at any given 
instant. This particular configuration would be equivalent to being on 
the lines designated as "B" on Fig. 3.13 through Fig. 3.15. 

These principal mode shapes can be roughly divided into 
acoustically and panel dominant types. The acoustically dominant modes 
are characterized by mode shapes and natural frequencies fairly close 
to the classical acoustic cases for a hard room. Figure 3.4 is an 
example of this type of mode. Figure 3.6 shows a panel dominant mode 
shape which has a natural frequency closely related to the "in vacuo" 
response for the panel and a highly disturbed acoustical mode shape. 
For each panel dominant mode there are two acoustical dominants. 

Further inspection of the mode shapes show that the nodal 
lines (zero pressure lines) have moved towards the flexible panel 
for the acoustically dominant modes and away from the panel for the 


panel dominant mode. This can be explained in terms of "mass" and 
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Figure 3.4 FIRST MODE - SOUND TRANSMISSION 
BETWEEN SYMMETRIC ROOMS 
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Figure 3.5 SECOND MODE - SOUND TRANSMISSION 
BETWEEN SYMMETRIC ROOMS 
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Figure 3.6 THIRD MODE - SOUND TRANSMISSION 
BETWEEN SYMMETRIC ROOMS 
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Figure 3.7 FOURTH MODE - SOUND TRANSMISSION 
BETWEEN SYMMETRIC ROOMS 
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Figure 3.8 FIFTH MODE - SOUND TRANSMISSION 
BETWEEN SYMMETRIC ROOMS 
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Figure 3.13 FULL COUPLING PANEL - RESONANT FREQUENCIES 
FOR SOUND TRANSMISSION 
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Figure 3.14 73.7 PER CENT COUPLING PANEL - RESONANT FREQUENCIES 
FOR SOUND TRANSMISSION 
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Figure 3.15 42.5 PER CENT COUPLING PANEL - RESONANT FREQUENCIES 


FOR SOUND TRANSMISSION 
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"stiffness" effects, as previously described. For any given en- 
closure the panel and the air in the adjacent room act as a "limp" 
(i.e. unstiffened) mass on the end of the enclosure to give the 
lower acoustically dominant mode. The second acoustically dominant 
mode occurs when the adjacent room reaches a frequency where it acts 
as a stiffener for the "limp" plate. Panel dominant modes occur at 
frequencies where the mass or stiffness effects of the rooms re- 
spectively decrease or increase the "in vacuo" panel frequency. 

It is noteworthy that Fig. 3.5 shows the paradoxical case 
of sound transmission with no panel motion. Though such a case can 
not exist with damping present, a natural frequency closely related 
to the one for this configuration was measured. 

Figure 3.9 through Fig. 3.15 relate the effects of changing 
receiving room length on the natural frequencies and principal mode 
Shapes for any given panel. Figure 3.9 through Fig. 3.12 were scaled 
with respect to the maximum pressure and the individual room lengths 
to give a basis for comparison of theoretical and experimental results. 
It is to be noted that good agreement occurs for all cases except 
where there is little or no panel motion, i.e. transverse modes and 
symmetric acoustical modes. These adverse configurations are highly 
dependent on damping which is not considered in this thesis. To the 
right of the scaled pressure diagrams are scaled displacement diagrams 
of the plan view of the coupling panel. 

The principal mode shapes within the coupled rooms can 
basically be described in terms of strain energy or frequency; the 


higher the strain energy stored, the higher the natural frequency 
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to maintain the energy level. Acoustically, the lowest strain 
energy state is the uniform pressure mode in which the room's 
overall pressure raises and lowers "en mass" while the second lowest 
is a single noded standing wave in the direction of the maximum room 
dimension. The panel is at its lowest potential state when unflexed, 
while the next higher state is a set of half-sine waves in the two 
planar directions, as previously described. When exciting a coupled 
system with any given frequency, the individual components will react 
more violently if the exciting frequency is near one of the components' 
natural frequencies. The resulting overall natural frequency wil] 
therefore be a modification of the components' natural frequencies. 

The paragraphs that follow describe each of the lowest five 
overall natural frequencies as a function of receiving room length 
starting with a short receiving room. In both the acoustical and 
flexural sense "mass" effect lowers a component's natural frequency 
while "stiffness" effect increases it. Also, included is some insight 
as to the effect of changing the panel area and stiffness on the over- 
all natural frequencies. The vertical lines "A", "B" and "C" in 
Fig. 3.13 through 3.15 will serve to locate zone points during dis- 
cussion. 

For the lowest overall natural frequency near point A, the 
receiving room has an almost zero-frequency uniform pressure which 
acts as a mass in tandem with the mass effect of the panel lowering 
the component frequency of the source enclosure. As the configuration 
moves towards "B" the acoustical stiffness of the receiving room drops 


to equal the stiffness of the source room; the wall still acts as a 
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limp mass, lowering the component frequency of either room. Moving 
towards "C" the receiving room develops a longitudinal single node 
standing wave, while the source room goes to a uniform pressure state. 
In this region the frequency is close to but above the component 
frequency for the receiving room. As the panel frequency increases 
and the panel forcing area decreases, the plate stiffness tends to 
raise the acoustic frequency of the system to the point where the 
source room appears to be acoustically hard and the frequency is that 
classically calculated (as in Fig. 3.15). It is to be noted that when 
the source room is shorter than the receiving room, the stiffness of 
the panel increases the frequency of the total configuration over that 
which would normally be calculated for the receiving room. 

The second mode is extremely dependent on the frequency of 
the coupling plate for unsymmetric rooms (points A and C). For the 
first two panel widths the natural frequency is close to the frequency 
for an "in vacuo" panel with a slight increase due to the acoustic 
stiffness of both rooms. When the rooms are symmetric, sound trans- 
mission with no wall motion takes place (totally independent of the 
plate frequency) at the natural frequency of the hard boundaried source 
or receiving room. For other configurations, acoustical mass effects 
of either room tend to lower the frequency of the panel. 

Beyond this point the overall natural frequencies and princi- 
pal mode shapes become disorganized. It is best to discuss Fig. 3.13 
and Fig. 3.14 and then Fig. 3.15 separately. 

The third principal mode for the two larger panel config- 


urations is governed by a single node longitudinal acoustic wave in 
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the receiving room in the vicinity of "A" and by the fundamental panel 
frequency between "B" and "C". The fourth and fifth mode shapes have 
transverse wave patterns in both source and receiving rooms regardless 
of receiving room length. In the fourth mode the wave patterns at 
any point near the coupling panel of one room is 180° out-of-phase 
with a point immediately adjacent to it in the adjacent room, causing 
a net force to be incident on any particular section of the panel. 
Since the fundamental frequency of the plate is below the overall 
natural frequency, the intercoupled plate-acoustic system has a 
natural frequency lower than that which would be calculated for simi- 
larly dimensioned "hard" rooms. The fifth mode has the adjacent wave 
patterns in phase so that the panel has no effect. 

For the more uncoupled configuration, Fig. 3.15, the second 
and third modes are degenerate in the region of "A". The frequency of 
the coupling panel is sufficiently above the frequency of a transverse 
acoustic wave that the previously independent frequencies associated 
with the "out-of-phase" and "in-phase" acoustic mode shapes are the 
same. This would indicate that the rooms are effectively hard. These 
mode shapes are similar to the fourth and fifth mode shapes for the 
wider panel configurations. In the region between "A" and "B" the 
second mode becomes dependent on the longitudinal acoustic mode in 
the receiving room and from "B" to "C" on the fundamental longi- 
tudinal acoustic mode in the source room. In the region of "B" the 
third and fourth modes are degenerate and remain as such for further 
lengthening of the receiving room. The fourth acoustic mode is de- 


pendent on the "in vacuo" panel frequency near "A" but becomes 
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72 
dependent on the fundamental longitudinal acoustic mode of the 
receiving room between "A" and "B". Beyond this, the acoustic and 
panel modes become coupled as the fundamental frequency of the panel 
is arrived at. 

Cases of "weak" intercoupling produce frequency curves that 
are almost discontinuous at points of curvature. These lines are 
almost co-incident with those calculated by "hard" room and "in vacuo" 
panel theories. One noteworthy exception is the longitudinal standing 
wave in the receiving room. The panel stiffness has some effect in 
increasing the lower frequency results and the mass effects of the 
panel decreasing the higher frequencies. 

In Fig. 3.13 through Fig. 3.15 the vertical distance between 
the pinned-pinned and pinned-fixed lines for any given mode is an indi- 
cation of the mode's dependency on the panel's boundary conditions. 
This exemplified itself when measurements were made on the experiment 
apparatus. Resonances associated with narrow distances between pinned- 
pinned and pinned-fixed boundary conditions were pronounced over a 
small band width whereas resonances associated with wider distances 
between the boundary conditions had a much wider "skirt". An adverse 
effect of the different "skirt" widths was the disappearance of 
certain less prominant resonances. Also, a shift occurred in the peaks 
of two closely adjacent resonances because of the phase plane addition 
of the respective skirt zones. Examples of this can be seen in 
Fig. 3.13 and Fig. 3.14 marked with "*". 

In Fig. 3.13, a 'W' denotes a resonance situation of the upper 


"rigid" panel of the receiving room and the coupling wall. It is given 
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Figure 3.16 42.5 PER CENT ORIFICE COUPLING - RESONANT 
FREQUENCIES FOR SOUND TRANSMISSION 
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Figure 3.17 4.2 PER CENT ORIFICE COUPLING - RESONANT 
FREQUENCIES FOR SOUND TRANSMISSION 
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because it appears as a resonance peak for the given configuration 


but is not sufficiently large in any of the other configurations. 


3.3 Sound Transmission Between Orifice Coupled Rooms 


Some interest was shown by Harris and Feshback [34] in 
the coupling of two rooms by a window. Considering the similarity in 
rooms studied, a slight diversion was taken to show the natural fre- 
quencies for a room with orifice coupling in a panel that represents 
43.7 per cent of the coupling wall. A second set of frequencies were 
calculated for an orifice coupling which represented 4.25 per cent of 
the coupling wall. These are shown in Fig. 3.16 and Fig. 3.17, re- 
spectively. 

For the large coupling window there is no extensive change of 
modal patterns from any given classical principal mode shape due to 
the excellent acoustic flow between the rooms. The constriction causes 
local pressure build-up in the corners adjacent to the coupling wall in 
longitudinal modes (shown in Fig. 3.18), which is common for both widths 
of constriction. The narrow window is effectively a line source of 
sound in the shorter room, whether it be the source or receiving room 
(shown in Fig. 3.18 (b)). 

From Fig. 3.16 it is reasonably easy to see that the coupling 
window has very little effect on the principal frequency. The natural 
frequencies for a similar room with no window are superimposed with a 
dashed line and are named. The deviation is due to the slight increase 
in kinetic energy as the flow passes through the common boundary. It 


is to be noted that the computational error is small, as indicated by 


re oe 
| 1 ht ae aor 


noi tequel tno 024 aitp-i ihe 


_. enottevwupt tio. * io wm 


nm 
ie 0 as 
i oe, es ae 
at [SE] Xosdri2s4 bos ef 18 wore 2: 
a 


: we 
live a 
t yvttesftmte ns mtn) 0 a 


7 ; 

= th 
-91t Sevwisn sit wode oF nodbt zon alii tne 
2ingesyqsy tedt fonsq 5 nt orth feel sof Tho Agr 
| 


oiew 29fonsupsyt to tae baoose A BW pat quay a 


os 


to tnea 19q @S.A betasesxge dorniw ent auo3! 82th 


-97 ,VI.€ .pFa bas Of. ert ah, nworlz-s om ste 


a 


=, 


to sprsio svrenetxe on 27 syst wobniw ont Tquon apne 


Ot sb. sGéde sbom [sqtantta [sole2sfs a} ye ct aniast 
— - 
282652 norsofyienoo sai .emooyr ont neswisd Wott ped 28 tn ety 


ni Piew pnaftlfquos sat of Fneasphs evsnybo Sad at qu-bi ted “g1u229" 
de : 


edtbiw rgod yot nomnoo ef dofdw .(87-8 .ofF AF nworta) a 


i ‘ ' . — Py ra 
fo soyuoe antl 5 ylavitostts 2t wobntw worsen ont 0 nttotT: 
P . a 7 
inno pnftyfeoey yo sovude edt sd 3F yerdtenw (Oey =p . si 2 ft 


a] 

® 

- 

— 

¥ ~ | 


paflavoo sat tent saz. ot yess yidanoesey 2f 3 art ‘etea 
; 


an 


yusen siT .yonsupeyt Teqronfiag odd no soetts SISSTT Yoav 
is A " aia 
itiw bezoqmrysque sys wobntw on ditw mooyr wei tmt2 he 
sesavont Snptla ort o3 gub-2F nottsiveb oat  .beasn off sant 
Si .yysbaued nomnoo sit Apuoyds zoc25q wot? ons 26 one 


(4 betsotbn’ 26 .[leme ef yore Tenorisduqmoo snd jens t 


ii 


the agreement of slope and only slight reversal of lines for the 
(3, 0, 0) mode. 

For the narrow passage, the natural frequencies and mode 
Shapes vary considerably from classical geometrically based calcu- 
lations. Again, the frequencies for the lowest mode are close to 
those for the (1, 0, 0) mode for a similarly dimensioned room with no 
coupling wall. Larger deviations occur due to greater kinetic energy 
through the orifice (much similar to two coupled Helmholtz resonators). 

For higher modes, there are indications of independence in 
room characteristics due to a smaller coupling surface. Modes three 
and four are both transverse wave patterns with the wave "in-phase" 
and "out-of-phase", respectively, much similar to the panel-coupled 
rooms. The greater frequency associated with mode four is due to a 
Slight pressure gradient that has its peak at the orifice and declines 
towards the end walls. 

Modes two and six also display this partially uncoupled 
characteristic. When the receiving room is short the receiving room 
has a transverse wave pattern, while the source room has a longitudinal 
standing wave. As the receiving room lengthens, the transverse wave 
in the receiving room is gradually replaced by a longitudinal wave in 
mode two. For mode six the transverse wave does not disperse but a 
longitudinal wave develops in the receiving room. 

Mode five has a longitudinal wave pattern much similar to the 
(3, 0, 0) mode of the equivalently dimensioned uncoupled configuration. 
As in mode one, the frequencies are lower because of the Helmholtz 


resonator effect. 
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CHAPTER IV 
DISCUSSION AND RESULTS 


There are several interesting effects which occur when 


the natural frequencies of the coupling plate are within the range 


of the natural acoustic frequencies. 


1) 


2) 


Pressure levels in the receiving room can 

be higher than those in the source room. 
Acoustically dominant modes occur at fre- 
quencies near those classically calculated 
with only slight deformation of the acoustic 
mode shapes. In this situation the exciting 
frequency is sufficiently removed from the 
natural frequencies of the panel so that the 
panel acts as a "passive" spring (stiffness) 
or "limp" mass on an acoustical system. 

Plate dominant modes have natural frequencies 
that vary as much as 8 per cent from the 

"in vacuo" classical panel frequencies, though 
changes in plate mode shapes are not appreciable. 
The mass and stiffness effects of the adjacent 
acoustical systems have considerable effect on 
the overall frequency, but subsequently suffer 
considerable deformation in acoustical mode 


shapes. 
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In the orifice coupled enclosures the effect of the 
aperture is not significant until the opening is less than 10 per cent 
of the common wall area. In this range the Helmholtz resonator effect 
causes the natural frequency to drop considerably, while the rooms 
become acoustically uncoupled to some extent. As the opening becomes 
smaller the straight isobars of a classical hard boundaried room be- 
come distorted into curves because the orifice becomes a stronger "line" 
source radiating into the room that is the least acoustically exicted. 

For both types of coupling, the insertion of an intermediate 
surface, be it a hard boundary with an air slit in it or a coupling 
panel, causes pressure build ups at these surfaces. This is due to 
the increase in potential energy and the decrease of kinetic energy 
because of reduced transverse velocities at a wall. 

The principal mode shapes produced give an insight of the 
acoustic isobars and the deflected form of the panel for a given resonant 
condition. With the mathematical characteristics of linear combination, 
visualization as well as calculation of composite sound fields is readily 


possible. 


4.2 Future Research 

There are some other possible extensions of the Kantorovich 
approach in plate-acoustic systems. One would be to study the random 
dynamic response of the aforementioned rooms and structure. Another 
would be to change coordinate systems. By developing a cylindrical 
acoustic element, a circular plate element and a cylindrical shell 


element (with the ©-direction being used for Fourier approximation) 
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and the required coupling matrices, the dynamic or resonant response 
of submerged closed cylindrical sections could be studied. The form 


of the approximating functions would be 


r" F (r,z) (A cos n © +B sin n 0) 


which would lend easily to the approximation of the Bessel Functions 
used in the solution of classical equations. Notably there would have 
to be care taken to circumvent the singularity that would occur at the 


origin (r = 0) of the element configuration. 


4.3 Results 

The finite elements produced give a good approximation of 
the principal mode shapes and natural frequencies for a coupled plate- 
acoustic system with two parallel surfaces and perpendicular sides. 
The agreement is best when the configuration is least dependent on 
internal hysteretic and viscous damping. Convergence of the approxi- 


4 or Vn,” for 


mate results to the exact solution is at a rate of I/n, 
all cases studied. The use of Kantorovich's approximation and Iron's 
symmetricization procedure produced a significant saving in computing 
time and storage. 

Such an approach is sufficiently general to allow study of 
many acoustic panel configurations with little alteration to the major 
portions of the computing program. If the program is to be used for 
various problems of varying geometry, the finite element approach is 


much more attractive than other analytical methods (i.e. finite 


difference). 
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uU,V,w - displacements in the x,y,z 
directions respectively 


A initial and final specific 
volumes 
el ae initial and final densities 


dx,dy,dz - dimensions of the infinitesimal 
element in the x,y,z directions 
respectively 


Figure Al INFINITESIMAL ACOUSTIC ELEMENT 
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APPENDIX A 
ACOUSTIC FORMULATION 


A.1 Equations of Motion for Acoustics 

A derivation of the equations is given by H. Lamb [8], 
therefore what is presented here will only represent a condensation 
of the above author's work. | 

Considering the element as shown in Fig. Al (with the 
included notation) and defining the dilatation A as the ratio of the 


specific volume increment to the original specific volume, we get 


A = (v - Vo) Z V9 
V = Vo (1 + A) (Al) 


Similarly, condensation s can be defined for densities 
thee (0 = Po) / Po 
Ben il tas) (A2) 
The density is the inverse of the specific volume, therefore 


T= Cl 3) (Fs) 


Considering only small dilatations and condensations, it can be 


shown that 


A = =s (A3) 


Looking more closely at the dilatation, it is defined by 
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Neglecting second order terms, it becomes 
where 


Poisson and Laplace consider the thermodynamics of sound to be almost 


adiabatic, which agreed well with experiments. Therefore, from (A2) 
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where y is the ratio of the specific heats. Again, assuming a small 


dilatation 
P = Py (1 + ys) (A5) 


By taking a force balance in each of the principal directions, we 


get 
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with the minus sign indicating the direction of acceleration 
relative to the pressure gradient. Equation (A6) can be written 


in more efficient notation as 


2 
- 06 
p a = grad p (A7) 
It is worthwhile to note that the velocities in the different 
coordinate directions within the acoustic element can be directly re- 


lated to the condensation s by 
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In acoustics the excess pressure (p - Po)» is used rather 
than the overall pressure p. Rearranging (A5) to get the excess 


pressure Pp), and making substitutions from (A3) and (A4) 


Das Baa, div u (A8) 


The excess pressure can be directly substituted into 


(A6) and upon combination with (A7) 


div (grad P,) Spey! P| “= 0 


The divergence of the gradient of a scalar is the 


Laplacian operator ve 


which is a standard form of the wave equation 
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and YP is the square of the wave velocity, which from now on is c, 


fe 
the velocity of sound 


sar iar (A10) 


A.2 Energy Equations for Acoustics 


The kinetic energy T for an infinitesimal element can be 


written as 


The potential energy U can be found by integrating the shaded area in 


Fig. A2 over the whole volume V 


U = SSS {Pg(V-Vp) - 1 (Pg - p) (v-vp) } dV 
Vv V Z Vv 
0 0 
Since the average volume change over the whole 7s zero, the first term 
disappears and by using (Al), (A3), (A8) and (A10), it can be written 


that 


U = pc® SIS (div u)® av 
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To this point no assumptions have been made as to the form of 
wave motion. Gladwell and Zimmerman [28] assumed a periodic sound wave 


and continued onto the energy equations. Craggs [15] studied the tran- 
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Sient response and, therefore, his formulation had to be more 
general but is, in some respects, easier to follow and will be used 
for the time being. 


A velocity potential » must be defined 
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where 2 is the value of the integral at time, ty. 
It can be derived from the previous equations without much 


difficulty that 


6u = grad ¢ 
ét 
divu=1 69 
ce ét 
P, = -9 oo 
ét (All) 


Note that from here on a superscript "+*" indicates differentiation 
with respect to time. The value of the potential lies in its ability 
to define the total character of a vector quantity u by a scalar, 


although at a loss in physical interpretation. 


A.3 Derivation of Equations of Motion from Functional Formulation 


It is now worthwhile to show that the first variation of 
the time integral of the Lagrangian, J, actually gives the correct 


equation. Making substitutions from (All) in the energy equation, 
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Applying Greene's theorem, the first term of 6J can be written as 
t 
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where S is the surface of the volume of interest. But 6d is zero at 
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must be zero. 
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These equations can be derived from (All) and represent, 
respectively, the wave equation and the natural boundary condition. 
It is only right that the "velocity potential" obey the wave equation 
as does the pressure p in (A9). The natural boundary condition can be 
translated either as the normal pressure gradient or the normal velocity | 


being zero at a rigid boundary. 


A.4 Matrix Formulation of the Functional 


Let the velocity potential within an element be written 
{o,(t)}, (A13) 


where {F(x)} is a column vector (of some length) of functions of x, 


the coordinates used to define the geometry, such that {o,(t) } (also 


e 
a column vector of the same length) is the value of » which varies with 
time at points called nodes. The right-hand side of (A13) may be 

differentiated and used in the desired equations containing st DUU Sie 
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neglecting the variables in round brackets. The energy equations 


for an element can now be written 
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which can be integrated to give the final form of the energy equations 
Pun T 
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The square of the acoustic density p is kept aside for consistency of 
units in later work where [S.J]. [PJ are the elemental acoustical 
stiffness and pressure matrices. These matrices can easily be defined 
by comparison of (Al4) and (A15). 
Substituting these equations in the functional 
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Taking the first variation of this equation and integrating the second 
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The first term is zero because §@ is zero at ty and to. This can be 


seen from the potential energy equation in (A1l4), knowing that 


60 = {F(x)}' {66}, = 0 at t, = O and t, = 0 


Since 6® is also arbitrary between ty and ty 
[S.J ese ag [PJ {onto = 0 (A16) 


The velocity potential lacks physical significance so rather than using 
the nodal vector to} (A1l6) can be differentiated with respect to time 
and rewritten in terms of the element nodal excess pressures {Pte by 


using (Al1) 
[SED Sh aL liner aes./0 (Al7) 


This last formula will be used in the finite element analysis 


of acoustic enclosures with rigid boundaries. 
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APPENDIX B 
MATHEMATICAL TABULATIONS AND OPERATIONS 


B.! Fourier Approximation and Integrals Over One-Half Period 


Let h(z) be a function that satisfies the wave equation 


and boundary conditions that sh = 0 at z=Oandz=d. If it is 
6Z 
piecewise continuous and finite within the interval, then 


where 
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The following integrals will be useful in future analysis 
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B.2 Hermitian Interpolation Polynomials 


i) Polynomials: O<e<] 


ii) Characteristics: 
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iii) Integrals: 
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B.4 Acoustic Rotational Matrix 
Since the element is orientated so as to simplify the surface 
integral, it must be rotated to a position where the nodal characteristics 
are compatible with adjacent elements. The symmetry of the [P] and [S] 
matrices can be retained if the invariance of the energy equations are 
used. The rotation to a new system (x', y') from the old system 


(x, y) iS a simple one: 
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To rotate the total pressure vector for the triangular 
element, three such rotational matrices must be placed down the 


diagonal of a large matrix, [0], whose remaining entries are zero 
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B.5 Coordinate Relation Matrix for Triangular Acoustic Elements 
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APPENDIX C 
MATRIX FORMULATIONS 


C.1 Decreasing the Order of a Matrix 


The solution to the eigenvalue problem for a single trianqular 
acoustic element gives ten eigenvectors in}. of which {nhig is of the 
least interest for this analysis (because of its high eigenvalue). These 
eigenvectors are orthogonal through the [S.] or [hal matrices of the 


equation 
2 x 
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The orthogonality of the eigenvectors allows any pressure 


{p,} within the element to be defined by 


By adding the constraint that the tenth eigenvector does 
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not contribute to the pressure distribution, a new pressure vector 


pee is arrived at such that 


{p,*}' [A] {n}y9 = 0 (C1) 


where [A] replaces the [S.J or [Pel 


By partitioning the vectors and matrix as follows 


p,* ny 
Po* Ny 
{pee = e {nhig 2 é 
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where Ay» Alo are vectors and Ano is a scalar. 
Let ip eet be a vector of the upper nine elements of ipeee and 


nohig be the upper portion of inhig: Then (C1) can be written as 
[ORE eR (nae (pee Dao tn tas 
rece onis “eltg ec 2! ‘o 


is T ig 


7 nia a 7 U ; a 
1 4 ; »@ gay 

; ‘Vee ay? 

soins swuezorg wan aot ra a“ | ot of oty 


80 f 


lavas 


af LA ; aj ent 
Th, i ei a 


=) 

: a ’ ver" bode 

dave 36 beviv: 
7 : 1 ] 


| 


| . isaac Co a yb Aan _ 
(f3) 0: “ot fA. (+93 ths et i = 
q Pinte M A “9 i P nora 


“aa ee wt enetan 7 t " 


ewollo? 25 xtetem bob 2102 ay ai 


j ' sell 7: Lan. 
oi a 7 i Deus 
aa il Shinui -|= [A] : Sve i 
: ; <0 0) ere 
} ~ ae ne se ome 5 ly an 


7 ~i : | ot oe pe ps “i 


-7st92 6 at os brs erode 918 of Ci 3A 9796 lw 
bre {* 3? to. 2dnamels anti naqqu: antt to. votoey & 9d cc om a © 


26 nedtiw Sd N62) (fa) nett rae ¥0 roe reg 


ee is i i: | ee oe ‘eer - 
Se ele atthe 7 
on? Is Moat silat j UpgA a 
; aie A aa Hots wy 4 ft " mm . 4 ; 
ihe _ us| coal 
7 ? . 
mn) a one ait iH le 


7 c 
a | 
~— an # q re 


109 


Therefore 


{Pyo*}! = {p,}) [A] / B 


where 


A. = [[A,,] inohig b [Ajo] {nyt] 


Again, using orthogonality 
{p,*}! [A*] {p,*} = constant 


where [A*] can also be either of the elemental matrices. The matrix 


can be partitioned 
ey i aero), Moe Py Mey peony 
Po ‘al i @ 2c ec 12 10 
. {Pigs} [A,1*1 {Pa}, + {Py Q*} [A,5*] {Py 9*} = constant 
Substituting from (C2) 
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It is now seen that the vector {poh can be looked on as a condensed 
pressure vector while the result in outer square brackets is a con- 
densed acoustic stiffness or mass matrix depending on the definition 


Of As). 


C.2 Applying Boundary Conditions 


The overall matrices when initially assembled have no 
boundary conditions. The system must be "tied down" before the eigen- 
value problem can be performed, because the unconstrained matrices are 
Singular. 

Several of the nodal characteristics are related to other 
noaalachanacrer istics y12e.) op =1-"Ctn 1 op, Or are zero. (hewlocter 
characteristics can be enties by simply Late the rows and columns 
which they (the zero characteristics) multiply with, because they 
supply nothing to the energy of the system. The zero characteristics 
are then deleted from the overall nodal characteristic vector {p,}s 
leaving a "Shrunken" energy matrix and nodal characteristic vector. 
The other defined nodal characteristic {p*} can be separated from the 
undefined ones 1p by simply shifting rows and columns of the overall 
matrices so that the calculations remain the same, i.e. 

" [Bi] 


Ct et tsp 


p* pe 


i = constant 


where the left hand side is the original nodal characteristic vector 


and original matrix and the right hand side is the re-arranged form. 
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A relationship exists between {p*} and {p,} 
oo =r Raat? 
By partitioning the matrix similarly to the way shown in 
Appendix C.1, a new overall energy matrix [B*] can be found in terms 
of the unknown nodal characters 


[B*] = [B,,] + [8,5] [R] + [R]' [B.,] 
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